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Besides the structure of interactions within networks, also the interactions between networks are of the 
outmost importance. We therefore study the outcome of the public goods game on two interdependent 
networks that are connected by means of a utility function, which determines how payoffs on both networks 
jointly influence the success of players in each individual network. We show that an unbiased coupling 
allows the spontaneous emergence of interdependent network reciprocity, which is capable to maintain 
healthy levels of public cooperation even in extremely adverse conditions. The mechanism, however, 
requires simultaneous formation of correlated cooperator clusters on both networks. If this does not emerge 
or if the coordination process is disturbed, network reciprocity fails, resulting in the total collapse of 
cooperation. Network interdependence can thus be exploited effectively to promote cooperation past the 
limits imposed by isolated networks, but only if the coordination between the interdependent networks is 
not disturbed. 



Not only are our social interactions limited and thus best described not by well- mixed models but rather by 
models entailing networks' it is also a fact that these networks are often interconnected and indeed very 
much interdependenf *" From the World economy to Google Circles, this interdependence cannot be 
denied and is all but absent, and in fact it is easily observable. Not necessarily so on the basis of actual links 
between the networks, which arguably could be the most hidden and difficult to identify feature of such systems, 
but much more so based on the impact actions in one network will have on the behavior in another network. 
Several examples are provided in', which demonstrate vividly that we live in a strongly interconnected World. 
More specifically, seminal works on interdependent networks have shown that even seemingly irrelevant changes 
in one network can have catastrophic and very much unexpected consequence in another network'''^. Since the 
evolution of cooperation in human societies also proceeds on such interdependent networks, it is thus of 
significant interest to determine to what extent the interdependence influences the outcome of evolutionary 
games. 

Preceding works concerning evolutionary games on individual networks and graphs, as reviewed comprehen- 
sively in' '" ", provide a stimulating starting point for explorations on interdependent networks'^ The discovery 
of network reciprocity due to Nowak and May'^, which demonstrated that the aggregation of cooperators into 
compact clusters can protect them against invading defectors even in the realm of the most challenging prisoner's 
dilemma game, launched a spree of research activity aimed at understanding the evolution of cooperation in 
structured populations. Methods of statistical physics in particular"", have proven very valuable for this task, as 
evidenced by the seminal studies of the evolution of cooperation on small-world'^'", scale-free"'^", coevolving^''^^, 
hierarchicaP-' and bipartite networks'''', to name only a few representative examples. Based on these studies it is 
now established that heterogeneous connections promote cooperation in all main social dUemmas^", and that this 
is a very robust evolutionary outcome^^, although not immune to the normalization of payoffs^*""^' and targeted 
removal of nodes'". In fact, heterogeneity in general, not just in terms of players having different degree within a 
network, has proven to be very effective in maintaining high levels of cooperation in the population-""''. 
Moreover, many coevolutionary rules" have been introduced that may generate such states spontaneously''""''*. 
Adding to this the discovered relevance of the population being structured rather than weU-mixed for the 
effectiveness of reward""'^", punishment^' selection pressure^'', and bargaining^^"^", it is beyond doubt that 
the added complexity of the interdependence between networks wUl lead to the discovery of fascinating new 
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insights concerning how and why we so often choose socially 
responsible actions over defection, despite the fact that large-scale 
social experiments suggest otherwise^'^. 

Here we wish to contribute to the subject by studying the evolution 
of public cooperation on two interdependent networks. In order to 
focus on the effects that are predominantly due to the interdepend- 
ence, we consider the square lattice as an elementary interaction 
network and an evolutionary game that is governed by group inter- 
actions, namely the spatial public goods game™. The usage of this 
setup is motivated by our desire to minimize other effects that could 
stem from network complexity or pairwise interactions. Neverthe- 
less, as we wiU argue when presenting the results and the newly 
discovered mechanism of interdependent network reciprocity, our 
findings should apply to many other social scenarios as well. In 
general, within the realm of evolutionary games, it is most natural 
to assume that the interconnectedness will not allow direct strategy 
exchanges between the networks, but rather affect the utility function 
of players which influences the strategy imitation probability within 
a network'^. While traditionally the utility is considered to be the 
payoff of a player as received from the games played with its inter- 
action partners, staging the game on interdependent networks opens 
up the possibility of alternative, possibly interdependent formula- 
tions of utility. A general example is when the payoff of a player is not 
just due to the games with its neighbors in the current network, but 
also due to its involvement in another network. Two individuals with 
strong ties but working in separate environments, or the relation 
between two companies of the same group but in different countries, 
can be characterized this way. We make use of these possibilities and 
propose a simple network- symmetric formulation of utility based on 
the payoff of player in the first network, as well as the average payoffs 
of its neighbors on both the original and the second network. While 
this formulation is not intended to model any specific scenario, it 
captures relevantly the essence of interdependent utilities. Notably, 
considering the payoffs of neighbors provides a natural coupling 
within a network, which can then be compared with the consequence 
of coupling between networks. Importantly, as the definition of util- 
ity is symmetric with respect to the two networks, there is no master- 
slave relation between them. As we will show, this setup leads us to 
the observation of interdependent network reciprocity. In particular, 
if related players on both networks aggregate into compact coopera- 
tive domains, then this form of reciprocity is much more powerful 
than traditional network reciprocity on an isolated network. On the 
other hand, even if some cooperators aggregate into a compact 



domain in one network, they ultimately fail to survive in the absence 
of similarly aggregated partners in the other network. Hence the term 
interdependent network reciprocity. 

Results 

We first present results for a = fl (see Methods for the description of 
the model and parameters), in which case the focus is mainly on the 
transition of importance from individual payoffs to average payoffs 
of nearest neighbors in the evaluation of utility. Naturally, as [i 
increases the interdependence between the two networks grows too 
(see Eq. 2), but it is impossible to attribute the reported effects 
unequivocally either to the growing influence of group payoffs or 
to the growing interdependence. As results presented in Fig. 1 show, 
simultaneously increasing both a and [i continuously lowers the 
minimally required multiplication factor r that is needed for coop- 
erator to survive. What is more, the span of r where cooperators and 
defectors coexist shrinks too, ultimately leading also to an ever- faster 
arrival to the pure C phase as r increases further. More precisely, the 
minimally required r/ G drops from 0.748 for a = /J = 0 to r/G = 0.33 
for a = P = 0.45, and correspondingly, the required r/G for complete 
cooperator dominance drops from r/G = 1.09 to just r/G = 0.35. This 
is certainly impressive and motivating for wanting to understand the 
mechanism behind the remarkable promotion of public cooperation. 
Even before further results are presented, however, it is clear that the 
transition of weight from individual payoffs to average payoffs of 
nearest neighbors in determining the utility of each player plays an 
important role. In particular, the transition magnifies the benefit of C 
- C coupling on both networks, yet this effect alone would be unable 
to improve the viability of cooperators to such an extent as reported 
in Fig. 1. 

It is therefore of interest to decouple the impact of a and [i by 
eliminating the a = /i restrain, as presented in Fig. 2, where contour 
plots encode the stationary fraction of cooperators pc for three dif- 
ferent values of r/G. Regardless of the multiplication factor, it can be 
observed that at any given value of P, increasing a will elevate pc- 
Note that increasing a simply elevates the relevance of the payoffs of 
nearest neighbors on the expense of individual payoffs on any given 
network, without modifying the level of interdependence between 
them. The impact of /i at any given a, however, is less clear-cut. Note 
that increasing P, however, also elevates the relevance of the payoffs 
of nearest neighbors on the expense of individual payoffs, but on 
separate networks, thus increasing the level of interdependence 
between them. In general, increasing /? at a fixed value of a also 
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Figure 1 | Shifting weiglit from individual payoffs to average payoffs of nearest neighbors in the evaluation of utility promotes the evolution of public 
cooperation. Depicted is the stationary fraction of cooperators pcon both lattices in dependence on the normalized multiplication factor r/ G for different 
a = fi values, as denoted in the figure legend. Due to the payoff-coupling within and between networks, there is a remarkable drop of critical r/G values 
that are required for cooperators to survive and dominate on both networks. Results were obtained by means of Monte Carlo (MC) simulations using 
2XLXL = 2X 300 X 300 system size and 2 X 10^^ MC steps. Averages were determined from the last 10" MC steps. 
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Figure 2 | Shifting weight fi-om individual payoffs to average payoffs of nearest neighbors promotes the evolution of cooperation, but it is important 
whether the shift entails the average payoffs of neighbors on the host network (a) or the average payoffs of neighbors on the other, interdependent 
network (P). Depicted are contour plots showing the stationary fraction of cooperators pc in dependence on both a and /? for r/G = 0.4 (left), 0.7 
(middle), and 0.8 (right). Increasing a at any given /S promotes the evolution of cooperation. The same is true vice versa, only that the positive effect 
depends much more on a. and r/G, thus highlighting the importance of an intact network reciprocity on each individual network for increasing ft to work. 
Details of Monte Carlo simulations are the same as described in Fig. 1. 



elevates pc levels, thus suggesting that an increase in the inter- 
dependence of the utility function also promotes public cooperation, 
similarly as the increase of a at any given Yet the efficiency of 
increasing ji depends much more on a than vice versa. For low values 
of a increasing ji has a much more marginal impact than if increasing 
P at high values of a. This difference in impact also grows as r/G 
decreases. Both these facts indicate that network reciprocity on each 
individual network is crucial for higher (i to take effect, as both lower 
values of a and r/G make it increasingly difficult for cooperators to 
form compact domains. 

Having then established the relevance of traditional network reci- 
procity on each individual network, it is next instructive to monitor 
the evolution of cooperation separately on the two networks. In order 
to do so, we first examine how the fraction of C - C and D-D pairs 
varies separately on network A and B for a characteristic combina- 
tion of game parameters, as presented in Fig. 3. Here P, and denote 
the fraction of the related pair configurations within and between 
networks, respectively. Accordingly, P,(C,C) is the probability to find 
a C - C pair in any given network, while Pg{C,C) is the probability to 



find a C - C link between the two networks. Besides the initial fall of 
cooperators before the eventual rise, which is the hallmark signature 
indicating the spontaneous emergence of spatial reciprocity, it is 
important to note just how synchronized the rise (fall) and fall (rise) 
of D - D (C - C) pairs on the two networks is. There is only a small 
quantitative difference inferable, but otherwise both cooperators and 
defectors on the two networks share completely the same fate. 

The question that now poses itself is just how important the syn- 
chrony of the evolutionary process on the two interdependent net- 
works is. As evidenced by the results presented in Fig. 4, the 
importance can hardly be overstated. If we prevent the spontaneous 
formation of cooperative domains in one of the networks by applying 
a permanent mixing of strategies, then this will negatively influence 
the evolution of cooperation in the other (undisturbed) network as 
well. On the latter, the cooperators will die out despite the fact that 
the coupling within the network is still functioning via nonzero a, 
and that thus the evolution of cooperation ought to be supported on 
this imdisturbed network. As can be observed, this intervention on one 
of the networks will not only prevent the emergence of synchronized 
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Figure 3 | Cooperators and defectors on the two interdependent networks share exactly the same fate during the course of evolution, i.e., as soon as the 
density of cooperator (defector) pairs starts decreasing (increasing) on one network, the same happens on the other network. In the legend P, and Pj 
mark the corresponding pair configuration probabilities within (internal) and between (external) networks, respectively. The evolution is completely 
synchronous, with quantitative differences emerging only due to the applied relatively small value of ft. Notably, for ft 0.5 the difference would be 
virtually none. Parameter values used were: r/G = 0.7, a = 0.15, and ft = 0.5, at L = 300. 
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Figure 4 | Mixing the distribution of players in one of the two networks disrupts the synchronous evolution of the two strategies, which has dire 
consequences for the fate of cooperation on the other network, despite the fact that the coupling of payoffs within the network is stUl significant. 

The evolution of cooperation in the latter, measured by the density of cooperators, is denoted dashed red. For comparison, solid green line shows the time 
dependence in the same network while the strategy evolution is not disturbed in the other network. This experiment highlights that while traditional 
network reciprocity alone can be inefficient, the chances of cooperators improve dramatically if interdependent network reciprocity is allowed 
to take effect. The applied system size was 2 X 200 X 200 at r/G = 0.7, a = 0.15, and ft = 0.5. 



evolution on the two networks, but will also result in fuU defection in 
the undisturbed network in spite of the fact that the two networks are 
not physically connected (the corresponding evolution in the unper- 
turbed network is denoted by "disturbed" in Fig. 4). For comparison, 
we also plot the evolution in the same network when the emergence of 
parallel evolution was not disturbed (denoted by "free" in Fig. 4). As 
the "fi-ee" line shows, cooperators can eventually spread widely even at 
parameter values where traditional single-network reciprocity would 
fail completely. 

These results lead us to the coining of the term "interdependent 
network reciprocity", simply because if network reciprocity fails on 



one network, it wiU inevitably fail also on the other network. If it 
remains intact on both, however, and if the interdependence via the 
utility function is sufficiently strong to initialize spontaneous coor- 
dination between the two networks, then interdependent network 
reciprocity will strongly outperform traditional single-network reci- 
procity, in fact leading to much higher levels of public cooperation, as 
evidenced in Figs. 1 and 2. The mechanism of interdependent net- 
work reciprocity can be demonstrated most beautifully by means of 
characteristic snapshots of the two interdependent networks as 
obtained from a prepared initial state. Figure 5 shows that only the 
one circular cooperative domain is able to survive the initial 




Figure 5 | Direct observation of interdependent network reciprocity. From left to right, panels (a)-(d) depict the evolution of cooperation in one 
network, while panels (e)-(h) show the same for the other network. Starting from a prepared initial state, only those circular cooperative domains that are 
initially present on both networks, and which can therefore make immediate use of interdependent network reciprocity, survive and eventually spread 
across both networks. Cooperators that are initially distributed uniformly at random, as well as clustered cooperators present on one but not the other 
network, surrender under the evolutionary pressure from defectors already at the early stages of evolution. The applied system size was 2 X 200 X 200 at 
r/G = 0.66, o; = 0.15, and ^ = 0.5. Snapshots were taken at t = 0, 100, 500 and 10000 MC steps. 
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Figure 6 | Interdependent network reciprocity works also one more than two symmetrically interdependent networks. Like for two networks, for three 
interdependent networks shifting weight from individual payoffs to average payoffs of nearest neighbors also promotes the evolution of cooperation, but 
it is important whether the shift entails the average payoffs of neighbors on the host network (a) or the average payoffs of neighbors on the other, 
interdependent network (/J) . For interdependent network reciprocity to work best, a threshold value of a ~ 0. 1 7 needs to be surpassed. The applied system 
size was 3 X I X L = 3 X 400 X 400 and 5 X 10^^ MC steps were used, while the value of r/G = 0.7 was the same as used in the middle panel of Fig. 2, which 
thus also serves as a comparison reference. 



onslaught of defectors, namely the one that is initially present on both 
networks, but not the other two which are initially present on only one 
network. Notably, other cooperators that are initially distributed uni- 
formly at random, and which could potentially be saved by traditional 
single-network reciprocity, die out fast too. Thus, only those coopera- 
tors that can make immediate use of interdependent network recipro- 
city survive and eventually spread across both networks. 

Lastly it is of interest to elaborate on the generality of the identified 
interdependent network reciprocity. To do so, we consider an 
alternative setup entailing three symmetrically interdependent net- 
works with the governing utility function being 

U, = {l-a-2li)P, + a(P,) -h li{P,,) + li{P,.,). (1) 

Here x' and x" denote the related players of player x in the corres- 
ponding networks. In this way every network is connected with two 
other networks in a symmetric way, without any master-slave pre- 
ferences. As in the originally studied model, strategy imitation is 
allowed only between players residing on the same network, but 
never between players from different networks. As results presented 
in Fig. 6 show, the impact of different a and fi values is qualitatively 
the same as reported in Fig. 2 for the two-network setup. In particu- 
lar, it can be observed that at any given value of /J, increasing a will 
elevate pc> while increasing ji at a fixed value of a also elevates pc 
levels, but it requires a threshold value of a to be surpassed for the 
interdependent network reciprocity to become really effective. 
Altogether, these results indicate clearly that the spontaneous emer- 
gence of coordination and subsequent synchronization of the evolu- 
tionary process is possible also if more than two networks are 
interdependent. Further on the generality of interdependent network 
reciprocity, we remind that the outcome of public goods games is in 
general independent of the interaction topology''", and therefore a 
similar impact can be expected also if the topology of individual 
networks is different from the square lattice. Moreover, since social 
dilemma games governed by pairwise interactions were recently also 
found to be susceptible to the "wisdom of group" effect""', we argue 
that interdependent network reciprocity should be observable 
straightforwardly by them as well. More challenging setups, which 
would be of interest to study more accurately in the future, would 
entail determining to what degree the network topology of each 
individual network forming the interdependent system is allowed 



to differ before for interdependent network reciprocity to remain 
intact. 

Discussion 

Summarizing, we have studied the evolution of cooperation in the 
public goods game on interdependent networks that are subject to 
interconnectedness by means of a network- symmetric definition of 
utility. Strategy imitation has been allowed only between players res- 
iding on the same network, but not between players on different net- 
works. We have shown first that, in general, increasing the relevance of 
the average payoff of nearest neighbors on the expense of individual 
payoffs in the evaluation of utility increases the survivability of coop- 
erators. More importantly, however, we have shown that while doing 
this for the coupling between individual and group payoffs on a given 
network simply reinforces network reciprocity, applying the same 
coupling of payoffs between networks has a much more interesting 
and in fact imexpected outcome. In particular, we have identified the 
spontaneous emergence of the so-called interdependent network reci- 
procity. The later impUes that if cooperative players on network A 
aggregate into compact domains while the corresponding players on 
network B fail to do so (or vice versa due to the network-symmetric 
definition of utility), the mechanism of network reciprocity will ulti- 
mately fail on both networks. On the other hand, if players on network 
A aggregate and the corresponding players on network B do too, then 
the evolution of cooperation wiU be promoted much stronger than 
based on traditional single-network reciprocity alone. At this point it 
is important to note that players between the networks are not phys- 
ically connected, nor are they allowed to imitate strategy from one 
another. The coordination of clustering of cooperators on both net- 
works is a spontaneous process that is thus brought about solely by the 
interdependence of utility. We have shown the mechanism of inter- 
dependent network reciprocity to be robust against the number of 
interdependent networks, and we have provided arguments in favor 
of robustness against variations of the structure of each individual 
network, as well against the type of the governing evolutionary game 
played. Overall, our research indicates that the interdependence of 
interaction networks can be exploited effectively to promote pubUc 
cooperation past the limits imposed by isolated networks, although 
this is subject to warranting undisturbed coordination of the distri- 
bution of strategies on the two networks, which in turn requires a 
sufficiendy strong interdependence of utility. 
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Methods 

The public goods game on both networks is staged on a L X L square lattice with 
periodic boundary conditions, where players are arranged into overlapping groups of 
size G — 5. Every player is thus surrounded by its ~ G — 1 nearest neighbors and is a 
member in^ — G different groups. Initially each player on site x in network A and on 
site x' in network B is designated either as a cooperator C or defector D with equal 
probability. The accumulation of payoffs and Py on both networks follows the 
same standard procedure. Namely, in each group cooperators contribute 1 to the 
public good while defectors contribute nothing. The sum of contributions is subse- 
quently multiplied by the factor r > 1, reflecting the synergetic effects of cooperation, 
and the resulting amount is equally shared amongst the G group members. In each 
group the payoff obtained is Pf on network A and P^, on network B, while the total 

amount received in all the groups is thus P^ — X] ^x! g^x' ■ 

While the two networks are not physically connected, interdependence is intro- 
duced via the utility function 

U, = {l-^-P)P, + a{P,)+P{P,,), (2) 

where (P^) and {P^') are the average payoffs of all four neighbors of players x and x' on 
their host lattices, respectively. It is important to note that this definition is network- 
symmetric. The evaluation of U^' is identical if the indexes ( ' ) for the two networks are 
switched, and thus there is no master- slave relation between them. The later was 
considered previously in^^, but it constitutes and entirely different setup. In accord- 
ance with the proposed symmetric definition of utility as given by Eq. 2, increasing a 
and/or ^, within the limitations of the condition {a + ji) < 1, puts more emphasis on 
the payoffs of the neighbors, while at the same time decreasing the relevance of 
individual payoffs. As expected, larger values of a promote cooperation, but the 
consequence of increasing ^ is hardly foreseeable as there is no direct information 
exchange between the networks via strategy imitation. Setting ^ — 0 decouples the 
two networks and more specifically, at a ~ ^ — 0 we regain the traditional spatial 
public goods game on both square lattices, where, as reported in^", cooperators die out 
if the multiplication factor r decreases below a threshold, depending also on the 
uncertainty governing the strategy adoptions K (see Eq. 3 below). For K ~ 0.5, which 
we will use throughout this work without loss of generality, the critical value is r/G — 

0. 748. 

Following the determination of utilities of players according to Eq. 2, strategy 
imitation is possible only between nearest neighbors on any given lattice, but never 
between players residing on different networks. We emphasize that this restrain, in 
addition to the network- symmetric formulation of utility, are the two key consid- 
erations of the present setup. Accordingly, on network A player x can adopt the 
strategy 5^ of one of its randomly chosen nearest neighbors y with a probability 
determined by the Fermi function 

W(Sy-^Sx)^ FJ-^ , , T , (3) 

l+exp[{Ux-Uy)/K]' 

where the utility Uy of player)' is evaluated identically as for player x. The probability 
of strategy invasion from player to player x' on network B is determined likewise, 
only that utilities U'^ and U'^ are used. Simulations of the model were performed by 
means of a random sequential update, where each player on both networks had a 
chance to pass its strategy once on average during a Monte Carlo (MC) step. The 
linear system size was varied from L — 100 to 600 in order to avoid finite size effects, 
and the equilibration required up to 10^ MC steps. 
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